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Abstract 

Electromagnetic (EM) wave scattering by many parallel infinite cylin- 
ders is studied asymptotically as a — > 0, where a is the radius of the cylin- 
ders. It is assumed that the centers of the cylinders distributed so 

that jV(A) = ln(-)J N(x)dx[l + o(l)], where A/"(A) is the number of points 

a A 

x m = (x m i,x m 2) in an arbitrary open subset of the plane xOy, the axes of 
cylinders are parallel to z-axis. The function N(x) > is a given continuous 
function. An equation for the self-consistent (limiting) field is derived as 
a — > 0. The cylinders are assumed perfectly conducting. Formula for the 
effective refraction coefficient of the new medium, obtained by embedding 
many thin cylinders into a given region, is derived. The numerical results 
presented demonstrate the validity of the proposed approach and its efficiency 
for solving the many-body scattering problems, as well as the possibility to 
create media with negative refraction coefficients. 

Key words: EM wave scattering by many thin cylinders; asymptotic so- 
lution; refraction coefficient; effective medium theory; nanowires; computa- 
tional modeling 

1. Introduction 



Wave scattering by many thin cylinders (nanowires) is important because 
of its many applications in chemistry [9], [22], medicine [23], optics [4], nan- 
otechnology [5], etc. Analytical formulas for solving electromagnetic (EM) 
wave scattering problem by many thin cylinders were derived in [19], and the 
results from [19] are used in this paper. 
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There is a large literature on EM wave scattering by arrays of parallel 
cylinders (see, for example, [7], [8]). Our approach has the following novel 
features: 

- the cylinders are thin, that is, they have small radius a, ka <C 1, where 
k is the wavenumber of the medium outside of the cylinders; this allows one 
to obtain a rigorous asymptotic solution of the wave scattering problem by 
many thin cylinders; 

- the solution to the wave scattering problem is considered also in the 
limit a — > when the number M = M(a) of the cylinders tends to infinity 
at a suitable rate and the distance d between neighboring cylinders is much 
greater than a, but there can be many small cylinders on the wavelength, 
so that the multiple scattering effects are essential; these effects are taken 
into account rigorously; both analytical and numerical methods for solving 
wave scattering problem on these thin cylinders (nanowires) are proposed 
and tested numerically; 

- the theoretical results obtained is a basis for a method for creating 
a new medium with negative refraction coefficient n 2 (x); this new medium 
is obtained by embedding many small perfectly conducting cylinders into a 
given (initial) medium. 

This work continues the earlier investigations in [10]-[21], and the numer- 
ical modeling presented in [l]-[3]. 

Let D m , 1 < m < M, be a set of non-intersecting domains on a plane P, 
which is xOy plane. Let x m e D m ,x m = {x m i,x m 2), be a point inside D m , 
S m be the boundary of D m , and C m is the cylinder with the cross-section 
D m and the axis, parallel to z- axis, passing through x m . We assume that x m 
is the center of the disc D m if D m is a disc of radius a. 

Let us assume that the cylinders are perfect conductors and a = 0.5diamZ), 
Our "smallness" (thinness) assumption is 

fca<l, (1) 

where k is the wavenumber in the region exterior to the union of the cylinders. 

We assume that the thin cylinders are distributed according to the fol- 
lowing law: 

Af(A) = In - / N(x)dx[l + o(l)], a -> 0, (2) 

CL J 
A 

where A/"(A) = 1 is the number of the cylinders in an arbitrary open 

& m eA 
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subset of the plane P, N(x) > is a continuous function, which can be 
chosen as we wish. The points x m are distributed in an arbitrary large but 
fixed bounded domain in the plane P. We denote by Q the union of domains 
D m , by Q' its complement in P. The complement in R 3 of the union C of 
the cylinders C m we denote by C 

The EM wave scattering problem consists of finding the solution to Maxwell's 
equations 

V x E = iu/iH, (3) 
V x H = -iueE, (4) 

in C such that 

E t = on dC, (5) 

where dC is the union of the surfaces C m , E t is the tangential component of 
E, \i and e are constants in C, u is the frequency, k 2 = uj 2 eji. Denote by 
Uq = so k 2 = uo 2 n 2 Q. 

Let us look for the solution to problem (3)- (5) of the form 

E(x) = E (x) + v(x), x = (x 1 ,x 2 , x 3 ) = (x, y, z) = (x, z), (6) 

where E (x) is the incident field, v(x) is scattered field satisfying the radiation 
condition 

^-&)=o(l), r = (x 2 1 + x 2 2 ) 1 / 2 , (7) 
and we assume that 

Eq(x) = e lKy+ik3Z e u k 2 + kl = k 2 , (8) 

where {ej}, j = 1, 2, 3, are the unit basis vectors along the Cartesian coordi- 
nate axes x,y,z. We consider EM waves with H 3 := H z = 0, i.e., -E-waves, 
or TiJ-waves, 

E = Y j E 3 e j , H = H iei + H 2 e 2 = ^^. (9) 

It is proved in [19] that the components of E can be expressed by the 
formulas: 
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E 3 = ^u x .e ik * z , j = 1, 2, E 3 = ue^ z , (10) 

K 

where u Xj := J^; the function u = u(x,y) solves the problem 

(A 2 + k 2 )u = in Q', k 2 := k 2 - fc 2 , (11) 

«lao = 0. (12) 

M = e ^ + w, (13) 

and it) satisfies the radiation condition (jTJ). 

One can check (see [19]) that the unique solution to (fTT]) - (fl3"]) is given by 
the formulas: 

E l = l -^u x e ik * z , E 2 = % ^u y ^\ E 3 = ue ik * z , (14) 

E x = ^u y e ik ^, H 2 = % ^u x ^\ H 3 = 0, (15) 

where u x := u y is defined similarly, and u = u(x) = u(x,y) solves scalar 
two-dimensional problem ( TTT]) -( |1~3|) . It is proven in [12] that such a problem 
has a unique solution. 

In [19] an asymptotic formula for this solution is derived as a — > 0. The 
results consist of the formulas for the solution to the scattering problem, of 
the equation for the effective field in the new medium obtained by embed- 
ding many thin perfectly conducting cylinders in the original homogeneous 
medium, characterized by the refraction coefficient n^, and of a formula for 
the refraction coefficient in the new medium. This formula shows that by 
choosing a suitable distribution of thin perfectly conducting cylinders, one 
can change the refraction coefficient, namely, one can make it smaller than 
Uq, and even negative. 

The paper is organized as follows. 

In Section 2 we derive a linear algebraic system (LAS) for finding some 
numbers that define the solution to problem ( Ilip - (IT5|) with M > 1, where M 
is number of cylinders. This is a new feature of our method: instead of looking 
for some unknown boundary functions (currents) we look for just numbers. 
This method is justified only if the cylinders are thin. We also derive an 
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integral equation for the effective (self-consistent) field in the medium with 
M = M(a) cylinders, M(a) — > oo as a — > 0. At the end of Section 2 these 
results are applied to the problem of creating a new medium with negative 
refraction coefficient by embedding many thin perfectly conducting cylinders 
into the original (initial) medium. 

In Section 3 the numerical results are presented. They demonstrate the 
validity and numerical efficiency of the proposed asymptotic method for solv- 
ing wave scattering problems. The relative error of the solution to the LAS, 
to which the wave scattering problem is reduced, is investigated; the optimal 
parameters M , a, and d that minimize the error of the asymptotic solution 
of the scattering problem are found. It is demonstrated numerically how the 
refraction coefficient of the new medium depends on the parameters M, a, 
and d. 

In Section 4 the conclusions are formulated. 

2. EM wave scattering by many thin cylinders 

In this Section, we derive LAS for the numbers u e (xj). These numbers 
determine the solution of the scattering problem by a rigorous asymptotic 
formula. Furthermore, we derive an integral equation for the limiting effective 
field, and obtain a simple explicit formula for the refraction coefficient n 2 of 
the new (limiting) medium. 

2.1. Asymptotic formulas for the effective field 

Let us assume that the domain D is a union of many small domains 

M 

D m , D = U D m . We assume for simplicity that D m is a circle of radius a 



centered at the point x m , and look for the solution to problem (11)-(13) of 
the form 



where S m is the boundary of D m , and dt is the element of the arclength of 



The distribution of the points Xjji 111 Sb bounded domain f2 on the plane 
P = xOy is given by formula ((2]). The incident field is Uq(x) := e lkX2 , and 



m=l 




(16) 
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g(x,t):=±Hg-\ K \x-t\). (17) 
The effective field acting on the Dj is defined by the formula 

u e = u^p = u(x) — J g(x,t)<7j(t)dt, \x — Xj\ > a, (18) 

Si 

or, equivalently, by the formula 

M 

u e (x)=u (x) + I g{x,t)a m {t)dt. (19) 

m=l,m^j Sm 

It is assumed that the distance d = d(a) between neighboring cylinders 
is much greater than a: 

d>a, lim-^- = 0. (20) 
d(a) 

Let us rewrite equation ffTBl) as 



M M . 

u = u + 9(x,x m )Q m + Y I [g{x,t) - g(x,x m )}a m {t)dt, (21) 

m=l m=l,m^jg m 

where 

Q m := J ° m {t)dt. (22) 

S rn 

It was proved in [19] that the second sum in (pPj) is negligible compared 
with the first one as a — > 0. The asymptotic formula for the numbers Q m is 
derived in [19]: 

Q m = ~ 2 ™f m) [l + o(l)], a^O. (23) 

a 

The new idea of our method consists of finding numbers Q m rather than un- 
known boundary functions o m (t) . This leads to a huge gain in the numerical 
efficiency of our method, and does not lead to the loss of its accuracy because 
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a is small. From formulas (21) and (23), one obtains the solution to problem 
(11)-(13) of the form, which is asymptotically, as a — > 0, exact: 



2tt m 

u(x) = u (x) - — p 9% x m )u e (x m ) + o{\). (24) 

The numbers u e (x m ), 1 < m < M, in (24) are not known. Below we derive 
LAS (25) and (31) for finding u e (x m ). The LAS (31) is of much lower order 
than the LAS (25), and can be interpreted as a collocation method for solving 
the integral equation (30) for the self-consistent (limiting) field in the new 
medium. The LAS (25), on the other hand, has a clear physical meaning. 

Setting x = Xj in (124"1) . neglecting o(l) term, and using the definition (1T91 
of the effective field, one gets a LAS for finding the numbers Q m : 



2tt m 

u e (xj) = u {xj) - y-j g{xj,x m )u e {x m ), 1 < j < M. (25) 

a m=l,mj^j 

This system can be easily solved numerically if the number M is not very 
large, say M < O(10 3 ). 

2.2. Integral equation for the limiting effective field 

If M is very large, M = M(a) — > oo, a — > 0, a linear integral equation 
for the limiting effective field in the new medium, obtained by embedding 
many thin perfectly conducting cylinders, is derived in [19]. 

Passing to the limit a — > in system (25) is done as in [16]. Consider a 
partition of the plane domain Q, in which the small discs D m are distributed, 
into a union of P small squares A p , of size b = b(a), a <^ b <^ d. For example, 
one may choose b = O^a 1 ^), d = (^(a 1 / 2 ), so that there are many discs D m 
in the square A p . We assume that squares A p and A q do not have common 
interior points if p ^ q. Let y p be the center of A p . One can also choose as 
y p any point X m 111 cL domain D m C A„. Since u e is a continuous function, 
one may approximate u e (x m ) by u e (y p ), provided that x m C A p . 

The error of this approximation is o(l) as a — > 0. Let us rewrite the sum 
in (25) as follows: 



2tt m p 1 

^ttt z2 g{x^x m )u e {x m ) = ln g(y q ,y P )u e (y p )^-j x > ( 26 ) 
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and use formula ([2]) in the form 



\ E l = N(y p )\A p \[l + o(l)]. (27) 



In , 



Here |A p | is the area of the square A p . 
From ff26l) and (127]) one obtains: 



27T P 



rE#j»^W^)= E 9(y q ,yp)N(y p )u e (y p )\A p \[l + o(l)}. (28) 



m a m^j p =l lP ^g 

The sum in the right-hand side of (128]) is the Riemannian sum for the integral 
^^9{y q ,y P )N{y p )u e {y p )\A p \ = J g(x,y)N(y)u( y y)dy, (29) 
where u(x) = limM e (x). Therefore, system fTSol) in the limit a — > yields 

a— i>0 

the following integral equation for the limiting effective (self- consistent) field 
u(x): 

w-^-tofaosmw. ,30) 

n 

One obtains a LAS for finding unknown quantities u(y q ), q = 1,2, P, 

see equation (|3T|) below, if one solves equation ( 130]) by a collocation method 
with piecewise- const ant basis functions. Convergence of this method to the 
unique solution of equation ( 130]) is proved in [13]. Existence and uniqueness of 
the solution to equation (13T?]) are proved as in [20] , where a three-dimensional 
analog of this equation was studied. 

The LAS (|3T]) is used for the numerical calculation of the limiting effective 
field and for a comparison of this solution with the solution to LAS (]25]) . 
whose order is much larger. The LAS is of the form: 



p 

u(y q ) = u (y q ) - 2tt 9(y q ,y P )N(y p )u(y p )\A p \, 1 < q < P. (31) 

Comparing the solution to ( 125]) with the solution to LAS ( 13T]) one finds the 
range of applicability of the asymptotic formula fliH]) for the effective field. 
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2.3. The refraction coefficient for the new medium 

Applying the operator A 2 + k 2 to equation ( 13"Uj) yields the following dif- 
ferential equation for u(x): 

A 2 u(x) + k 2 u{x) - 2nN(x)u{x) = 0. (32) 

This is a Shrodinger-type equation, and u(x) is the scattering solution cor- 
responding to the incident wave Uq = e lKy . 

If one assumes that N(x) = N is a constant, then it follows from (13"2"|) that 
the new (limiting) medium, obtained by embedding many perfectly conduct- 
ing circular cylinders, has new parameter k 2 n := n 2 — 2tcN. This means that 
k 2 = k 2 — k\ is replaced by k 2 = k 2 — 2nN. The quantity k 2 is not changed. 
One has k 2 = u 2 n 2 , k 2 = oo 2 nl. Consequently, n 2 /n1 = (k 2 — 2irN)/k 2 . 
Therefore, the new refraction coefficient n 2 is 

n 2 = n 2 Q (l - 2KNk~ 2 ). (33) 

Since the number N > is at our disposal, equation fl33|) shows that choos- 
ing suitable N one can create a medium with a smaller than n^, refraction 
coefficient, even with negative refraction coefficient n 2 . 

In practice one does not go to the limit a — )■ 0, but chooses a sufficiently 
small a. As a result, one obtains a medium with a refraction coefficient n 2 , 
which differs from (133]) . but the error tends to zero as a — > 0, and one has 
limn;? = n 2 . 

a-i-0 a 

3. Numerical results 

Some algorithms for computational modeling of the wave scattering by 
many small particles were developed in [1] for the acoustic wave scattering 
and generalized in [2], [3] for electromagnetic (EM) wave scattering. It was 
proved in these papers that the asymptotic solution of the many-body wave 
scattering problem, proposed in [10]- [16], [20], is computationally efficient and 
yields accurate numerical results. On the basis of this theory a recipe for 
creating a material with a desired refraction coefficient was formulated. This 
recipe was verified numerically in [2]. 

In this Section, numerical results are presented. These results demon- 
strate the efficiency of the asymptotical approach for solving the EM wave 
scattering problem in the media with many embedded perfectly conducting 
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cylinders of small radius a, and the possibility to create the medium with a 
negative refraction coefficient. 

The first portion of the numerical results demonstrates the approximation 
errors of the numerical solutions to the LAS ff25l) compared to the numerical 
solution of LAS f l3T]) . corresponding to the collocation method for solving of 
the integral equation f ]3"U]) for the limiting field. The rest of the numerical 
results demonstrate the possibility to create the media with negative refrac- 
tion coefficient, and yields optimal parameters M , a, and d, for creating such 
a coefficient. 

The following numerical problems are important from the practical point 
of view: 

- to determine the values of the parameters M, a, and d, that provide the 
solution to LAS (125]) with the desired accuracy, for example, with relative 
error of the order 10~ 3 — 10~ 4 ); 

- to investigate the convergence of the LAS fl3Tl) and to determine the op- 
timal values of the parameters M, a, and d, which provide such convergence; 

- to compare the solution to LAS (IZ51) and LAS f[3"T]) and to find the range 
of M, a, and d that provide accurate solutions; 

- to determine the values of M, a, and d, which yield negative refraction 
coefficient. 

3.1. The accuracy of the solution to LAS ( 125]) 

The numerical procedure for checking the accuracy to LAS ( I25p consists 
in calculations with various values of parameters M, a, and d at k — 1.41. 
First, we study the convergence of the solution depending on the number 
M of embedded into Q cylinders. The radius a of the cylinders is changed, 
the distance d between neighboring cylinders is kept in the range d > 10a. 
In Fig. 1, the relative errors of the solutions to (I25p are shown for the case 
when number M of cylinders grows from 25 to 3200. The values of \[M are 
depicted along the x axis. Because the exact solution to ( 125]) is not known, 
the relative error was calculated by formula 



RE 




(34) 



U 2 M 



instead of generally used 



RE 



(K 



mm |) 



(35) 
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where u 2 m and um are the solution to (I2"5j) with 2M and M cylinders respec- 
tively, is exact solution. 

The maximal error is observed for a/M = 5 and it is equal to 11.2%. The 
error diminishes when M grows and its value is 0.2% at M — 1600. The 
error practically does not depend on the radius a of cylinder for big M. The 
curves presented in Fig. 2 show that absolute error has similar character and 
its minimal value is 0.002 for considered a. 

The above numerical results are obtained for d > 10a. Computationally, 
the values of the error depend on this ratio. This can be used for minimization 
of the error by choosing different a and d. The numerical calculations show 
that there is an optimal value of the distance d between cylinders, which 
provide the minimal error if a is fixed. This optimal distance d depends on 
the number M of cylinders in D and varies in the range 14.5a — 40a. 




Figure 1: Relative error of solution to (25) versus M, k = 1.41 

Note that the errors of E components depend on the ratio of k and k 3 . The 
numerical results show that values of k (k 2 = k 2 — k 2 ) in a neighborhood of 1 
yield the minimal errors. This implies k 3 = 1.0 on account of k 2 = (1.41) 2 = 
2. 

3.2. The relative error for the solution to LAS (l3Tj) 

The collocation method [13] for solving LAS ( )3T|) . corresponding to the 
limiting equation (!30|) . is applied to check the accuracy of the numerical 
solution of LAS ( 13T|) . The relative error is defined in the previous subsection. 
In Fig. 3 the dependence of the relative error on the number P of the 
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Figure 2: Absolute error of solution to (25) versus M, k = 1.41 



collocation points is shown. When P is small, for example, P = 25, this 
error is large: it is equal to 30.2%, 27.4%, and 10.4% for a = 0.05, a = 0.01, 
and a = 0.001, respectively; the error is equal to 2% for a = 0.0001. In the 
considered range of P, the error depends on a. The value of ka does not 
exceed 0.0705 here. The smallest value of a = 0.0001 provides low error for 
all P. The values of absolute error are shown in Fig. 4. The maximal value 
of this error at y/~M = 35 does not exceed 0.005 and is achieved at a = 0.05. 




0.5 



5 10 15 20 25 30 V p 35 

Figure 3: Relative error of solution to (|3T|) versus P 

The values of the errors for the field u and its E components for the large 
P at a = 0.05 are shown in Table 1. The values of N(x) = N, with various 
constant values of N, are calculated by formula ([2]) and are shown in the last 
column of Table 1. 
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Abs.Err. 


Abs.Err. 


Abs.Err. 


Abs.Err. 


N 




of \u\ 


of \u\ 


of 


of \E 2 \ 


of \E 3 \ 




40 


0.76% 


0.0041 


0.0185 


0.0205 


0.0163 


61.53 


45 


0.54% 


0.0034 


0.0166 


0.0183 


0.0147 


77.89 


50 


0.36% 


0.0028 


0.0142 


0.0165 


0.0134 


96.16 


55 


0.22% 


0.0024 


0.0126 


0.0150 


0.0123 


116.36 


60 


0.10% 


0.0020 


0.0113 


0.0139 


0.0114 


138.47 



The above calculations were carried out at the values of that were 
determined by formula (jSJ) for the prescribed M and d. This implies the 
following formula 

"'A 

which agrees with formula (43) in [19] when N(x) = N = const, the dis- 
tribution of particles is assumed in the whole Q, and Q is the union of the 
non- intersecting domains A p . 

In formula (136]) the quantity Af(Q) is the total number of the embedded 
cylinders in Q, |A^| is the area of Q, k = 1.41. The calculations show that 
the value of N(x), calculated by formula fl36|) . can be varied so that it will 
provide the minimal error for the solution to LAS fl3T|) . In Fig. 5, the relative 
and absolute errors for the solution to LAS fl3~Tl) are shown in a neighborhood 
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of various values of N, calculated by formula ( 1361) . The first vertical line at 
the x axis corresponds to N = 6.49, calculated for a = 0.0001, and the 
second one corresponds to N = 7.86 for a = 0.0005. The minimal values 
of the errors are to the left of the values of N, calculated by fl36|) . For the 
considered parameters, the relative error decays from 0.72% to 0.31% for 
a = 0.0005, and it decays from 0.64% to 0.15% for a = 0.0001. The absolute 
error decays from 0.05 to 0.009 and from 0.025 to 0.005, respectively. 




2 4 6 S 10 N 12 



Figure 5: The absolute and relative errors of the field versus iV 




Figure 6: The absolute errors of E components versus TV 

The absolute error for the Ej, j = 1,2,3, components of the field is 
presented in Fig. 6. The errors for E\ and E2 are higher than that for 
E s because the components Ei and E 2 contain the derivative of function 
H^\kr) at the small values of kr, while E 3 does not contain the derivative. 
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The value of N plays a role of an additional parameter varying which 
one can decrease the error. The value of N can be changed by changing the 
distance d between neighboring cylinders while keeping fixed their number 
in the area where N is being changed. 

3.3. Comparison of solutions to LAS f l25|) and LAS fl3Tj) 

The accuracy of the asymptotic formula (24) was investigated by com- 
paring the solutions to LAS f )25|) and to LAS f )3T|) . The solution to fl3Tj) with 
P = 4900 collocation points is considered the benchmark solution, k = 1.41. 
The relative error of this solution does not exceed 1% at the considered values 
of a. This error is maximal at a = 0.05 and it decays if a decreases. 

In Fig. 7 the relative and absolute errors of the solution to LAS (25) are 
shown at various a. The maximal value of the relative error is observed at 
M = 25 and it is equal to 32.7%, 27.3%, and 16.6% at a = 0.01, a = 0.001, 
and a = 0.0001 respectively. This error for a/M = 35 is equal to 2.2%, 1.9%, 
and 2.4%. 

The absolute error for E% and E2 components is shown in Fig. 8. As in 
the preceding subsection (see Fig. 6) this error is higher than the error for 
E%. The largest error for E2 is equal to 1.2 for a = 0.01 when M = 25; the 
minimal value of the error for this M is obtained for E% at a = 0.0001, and is 
equal to 0.13. The minimal value of the error for E3 component is obtained 
when a/M = 35 and this error is 0.02. 




Figure 7: Absolute and relative errors of u versus number M of cylinders, 
k = 1.41 
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5 10 15 20 25 30 V M 35 

Figure 8: Absolute error of E components versus number M of cylinders, 
k = 1.41 

The above results are obtained in the case when the distance d between 
the cylinders is fixed. It turns out that the distance parameter d influences 
also the error of the solution to equation (24) when the number of the cylin- 
ders is fixed. The errors of the solution to equation (25) when M = 900 
for various values of d are shown in Fig. 9 and Fig. 10. The benchmark 
solution to LAS ( 130]) is the same as in the preceding example. There is an 
optimal value of d, which provides the minimal value of the error. The values 
of na, n = 1,2, 3, are shown along the x-axis. The minimal error equals 
to 0.57% and is obtained when a = 0.01 and na = 12.6; it is equal to 0.29% 
when a = 0.001 and na = 16.5, and it is equal to 0.23% when a = 0.0001 
and na = 21.6. The minimal values of absolute error when a = 0.001 and 
a = 0.0001 are shifted to the left in comparison to the minimal relative error. 

3.4. The refraction coefficient of the new medium 

One can conclude from formula (133]) that the value of the refraction co- 
efficient n 2 depends on the wave number k, on the parameter N, on a, on 
M, and on d. In Fig. 11 and Fig. 12 the dependence of n 2 on k is shown 
for two various values of d. The number M of the cylinders is equal to 225. 
The cylinders are placed equidistantly in a square 15 cylinders on the side of 
the square. The lengths I of square are equal to 0.1386m and 0.2857m in the 
Fig. 11 and Fig. 12 respectively. Consequently, the values of d are equal to 
0.0099m and 0.0204m. 
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Figure 9: Relative error of u versus distance d between cylinders, k — 1.41 
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Figure 10: Absolute error of u versus distance d between cylinders, k = 1.41 



The value of k has dimension L 1 , where L is length, the a and d have 
dimension L, and the values n 2 of the refraction coefficient are normalized 



to the value 11.1254 • 10 



-18sek 2 



This value is obtained by multiplying e 



8.85 x 10~ 12 F/m and /i = An x 10~ 7 H/m, taking into account the formula 
nl = e n , where F stands for farad, and H stands for henry, [H] = -y-, [•] 
stands for the dimension of a physical quantity, and T stands for time. 

At the smaller d = 0.0099m (see Fig. 11) the values of n 2 differ consid- 
erably from the refraction coefficient rig = 1 of initial media, because more 
cylinders are embedded per unit area. It is seen from Fig. 12 that the refrac- 
tion coefficient n 2 when d = 0.0204m is close to Uq. An increase of k forces 
n 2 to get closer to initial refraction coefficient n^. This is observed for all 
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considered values of a. The considered values of the M, a, and d yield the 
ratio 2 l ess than 0.05, so condition (I2"0"j) is satisfied. 
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Figure 11: Values of n 2 versus wave number k, M = 225, d = 0.0099m 
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Figure 12: Values of n 2 versus wave number k, M = 225, d = 0.0204m 

The numerical results presented in Fig. 13 demonstrate a possibility to 
create the medium with various refraction coefficients n 2 depending on the 
distance d between the cylinders, when a and k are fixed. The results are 
shown for a = 0.0001 and k = 20.0. At the small values of M the values of 
n 2 are changed considerably, and when M increases n 2 tends to the following 
values: n 2 = —0.45, n 2 = 0.06, n 2 = 0.35, and n 2 = 0.64 when d = 20a, 
d = 25a, d = 30a, and d = 40a, respectively. Note that at the considered 
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values of the parameters the relative error of the solution to LAS (25) does 
not exceed 2.34%, 1.69%, 1.18%, and 0.93% for d = 20a, d = 25a, d = 30a, 
and d = 40a at \/M = 30, and the relative error decays when M grows. 
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Figure 13: Values of n 2 



versus the number M of cylinders 



Consequently, one can change the refraction coefficient n 2 by changing k, 
a, M, and d. 

4. Conclusions 

Asymptotic solution is given for the problem of EM wave scattering by 
many perfectly conducting parallel cylinders of small radii a, ka « 1. An 
equation for the effective (self-consistent) field in the limiting medium is 
obtained when a — > and the distribution of the embedded cylinders is given 
by formula (j2J). The theory yields formula (1531 for the refraction coefficient 
of the new (limiting) medium obtained by embedding of these cylinders into 
the initial homogeneous medium. This formula shows how the distribution 
of the cylinders influences the refraction coefficient. 

The numerical results confirm the validity and efficiency of the asymptotic 
method for solving the above scattering problem. The optimal values of 
the parameters k, a, d, M, N, that minimize the error of the solution to the 
scattering problem, are found numerically. It is shown both theoretically 
and numerically that one can create negative refraction coefficients in the 
new medium. 
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